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ABSTRACT 

We present a self-consistent study of the stellar populations and the ionized gas in a 
sample of 10 HII galaxies with, at least, four measured electron temperatures and a 
precise determination of ionic abundances following the "direct method". We fitted 
the spectral energy distribution of the galaxies using the program STARLIGHT and 
Starburst99 libraries in order to quantify the contribution of the underlying stellar 
population to the equivalent width of H/3 (EW(H/3)), which amounts to about 10 % 
for most of the objects. We then studied the Wolf-Rayet stellar populations detected 
in seven of the galaxies. The presence of these populations and the EW(H/3) values, 
once corrected for the continuum contribution from underlying stars and UV dust 
absorption, indicate that the ionizing stellar populations were created following a 
continuous star formation episode of 10 Myr duration, hence WR stars may be present 
in all of objects even if they are not detected in some of them. 

The derived stellar features, the number of ionizing photons and the relative inten- 
sities of the most prominent emission lines were used as input parameters to compute 
tailored models with the photoionization code CLOUDY. Our models are able to ad- 
equately reproduce the thermal and ionization structure of these galaxies as deduced 
from their coUisionally excited emission lines. This indicates that ionic abundances 
can be derived following the " direct method" if the thermal structure of the ionized 
gas is well traced, hence no abundance discrepancy factors are implied for this kind 
of objects. Only the electron temperature of S"*" is overestimated by the models, with 
the corresponding underestimate of its abundance, pointing to the possible presence 
of outer shells of diffuse gas in these objects that have not been taken into account 
in our models. This kind of geometrical effects can affect the determination of the 
equivalent effective temperature of the ionizing cluster using calibrators which depend 
on low-excitation emission lines. 

Key words: galaxies: starburst - galaxies : stellar content - ISM: abundances - HII 
regions: abundances 



1 INTRODUCTION 

The understanding of the physical properties of gas, stars, 
and dust and their interrelations and evolution in differ- 
ent scenarios of star formation constitute one of the most 
intriguing open issues on astrophysics. This type of phe- 
nomenon is well studied at different spatial scales, since 
mEissive young star clusters are powerful sources of radi- 
ation, detectable even in the most distant galaxies. For in- 



stance, from a statistical point of view, large surveys of star- 
burst galaxies have been used to confirm, with a great confi- 
dence base, some of the observational relations between the 
physical properties of different types of galaxies at differ- 
ent cosmological epochs. This is the case of the luminosity- 
metallicity relation or its equivalent, the mass-metallicity 
relation, which establishes that in the local Universe the 
most massive galaxies have, in average, higher metallicities 
than the less massive ones (SDSS: Tremonti et al., 2004; 
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2dFGRS: Lamareille et al., 2004). This observational fact is 
a consequence of higher stellar yields in the more massive 
galaxies, which could find less troubles to retain the met- 
als ejected by stars and pushed out towards the interstellar 
and intergalactic medium by supernovae and stellar winds 
(Larson, 1974). 

Nevertheless, since this kind of studies are mostly based 
on the relation between the stellar mass and the oxygen gas- 
phase abundance, they present some drawbacks, involving 
several aspects in the determination of their properties, in- 
cluding the following: 1) The relation between the total mass 
of a galaxy and its stellar mass depends strongly on differ- 
ent environmental and morphological circumstances (Bell & 
de Jong, 2001). 2) The stellar mass is usually estimated by 
means of stellar synthesis population model fitting of the 
spectral energy distribution (SED). These models are of- 
ten affected by degeneracies between age, metallicity, and 
dust-extinction which are never well solved (Renzini & Buz- 
zoni, 1986). 3) The metallicity is only estimated in starburst 
galaxies through the oxygen gas phase abundance derived 
from the relative intensities of the most prominent emis- 
sion lines. This excludes from these studies the most qui- 
escent and early-type galaxies. 4) The spatial distribution 
of metallicity in massive, non-compact objects is often ne- 
glected, assigning a unique value to the whole distribution. 
However, it is known that the radial profile of the oxygen 
abundance in some close spiral galaxies can spread over more 
than an order of magnitude (Searle, 1971). 5) The strong- 
line parameters used to derive metallicities are often cali- 
brated employing sequences of models which do not follow 
the distribution of metallicity obtained with more confident 
observational methods like, for instance, the direct method 
(Perez-Montero & Diaz, 2005), which is based on the previ- 
ous determination of the electron temperature of the ionized 
gas by means of the fainter auroral emission lines. 

HII galaxies usually represent the low-end of the mass- 
metallicity relation in different surveys. In fact, the work of 
Lequeux et al. (1979), the first where the correlation between 
oxygen gas phase abundance and stellar mass was found, is 
based on this type of blue dwarf irregular galaxies. Their 
spectra are characterized by very bright emission lines while 
the total luminosity and color are dominated by the presence 
of intense bursts of recent star formation (Sargent & Searle, 
1970; French, 1980). They are dwarf and compact so they 
seldom show chemical inhomogeneities. The metallicity dis- 
tribution of HII galaxies has an average value (12-|-log(0/H) 
~ 8.0, Terlevich et al., 1991) sensibly lower than the rest of 
starburst galaxies and, in fact, the objects with the lowest 
metal content in the Local Universe belong to the subclass 
of HII galaxies (IZwl8, Searle & Sargent, 1972). 

Sloan Digital Sky Survey (SDSS) represents the largest 
sample of galaxies in the Local Universe and it has been 
explored for the study of the mass-metallicity relation 
(Tremonti et al., 2004). Unfortunately, this survey presents 
some drawbacks that prevent a precise determination of the 
metallicity in the emission lines-like objects: 1) Since the 
spectral coverage of the survey starts at 3800 A, the sample 
objects at redshift ~ have not any detection of the [On] 
A 3727 A line and hence no direct measurement of the total 
abundance of oxygen is possible. Only an estimate can be 
given in terms of the weaker auroral lines at AA 7319,7330 
A but with a high degree of uncertainty (Kniazev et al.. 



2003). 2) Only a fraction of the objects show in their spec- 
tra the weak auroral emission lines necessary to estimate 
the electron temperature and, hence, to calculate chemical 
abundances following the direct method (Izotov et al., 2006). 
3) The only auroral lines seen in the most part of this sub- 
sample is [Olll] A 4363 A and, therefore, only the electron 
temperature of the high excitation zone can be derived. The 
total ionization structure of the nebula is then often ob- 
tained based on this temperature through a very simplified 
assumption (Perez-Montero & Diaz, 2003). 4) For the rest of 
the objects with no direct electron temperature derivation, 
a strong-line method is used, so the determination of metal- 
licity is hugely sensitive to the chosen calibration (Kewley & 
Ellison, 2008). 5) The 3 arc sees aperture of the SDSS fiber 
precludes a good determination of the physical properties of 
the starburst region of these galaxies. 

Double-arm spectrographs are suitable tools to observe 
the entire spectral range between 3500 A and 1 micron in one 
single exposure at good spectral resolution the star forming 
regions of HII galaxies, avoiding the effects due to a different 
covering of the object in different spectral ranges. In the case 
of compact HII galaxies, their bursts of star formation can 
be covered entirely with a narrow slit. Assuming that differ- 
ent electron temperatures are measured in the optical spec- 
trum, this implies that different excitation regions inside the 
ionized gas can be characterized, leading to a more precise 
determination of the ionic chemical abundances. There are 
some reported discrepancies between abundance determina- 
tions obtained from bright coUisional lines and the fainter 
recombination ones (Peimbert et al., 2007), which do not 
depend on the determination of the electron temperature. 
Nevertheless, the existence of these abundance discrepancy 
factors (ADFs) is still not conclusively established. In fact, 
in the case of HII galaxies, no significant deviations have 
been found between the electron temperatures derived from 
coUisional lines and from the Balmer (Hagele et al., 2006) 
or Paschen discontinuities (Guseva et al., 2006). Besides, 
the high ADFs derived from optical recombination lines are 
sensibly lower than those calculated from mid-infra-red lines 
(Wu et al., 2008; Lebouteiller et al., 2008) or from planetary 
nebulae (Wilhams et al., 2008 ). 

In this work we model the stellar content and the abun- 
dances of the most representative ions of a sample of HII 
galaxies previously observed in the SDSS and re-observed 
with different double arm spectrographs in the spectral 
range between 3500 A and 1 micron (Hagele et al., 2006, 
hereafter H06 and Hagele et al., 2008, hereafter H08). The 
study of the stellar population in the star forming knots of 
these galaxies includes underlying and ionizing populations. 
We show that the correct characterization of the ionizing 
stellar population, including Wolf-Rayet stars, and an ap- 
propriate geometry of the gas, dust content and the abun- 
dances of the main elements are enough to reproduce with 
precision the intensities of the observed emission lines. We 
demonstrate too that it is possible to draw a consistent pic- 
ture of the thermal and ionization structure of the gas based 
on the measurement of as many electron temperatures as 
possible and, hence, to calculate ionic abundances with high 
precision using optical to far red coUisional emission lines. 

In the next section we describe the sample of objects 
that we have modeled. In section 3 we describe our results 
concerning the stellar content through model fitting of the 
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Table 1. List of the modeled objects, with some additional basic information. 



Object ID 


hereafter ID 




Publication 
(Telescope) 




redshift 


log L(Ha) 
(erg ■ s-l) 


12+log(0/H) 


Other names 


SDSS J002101. 03+005248.1 


J0021 




H06 (WHT) 




0.098 


42.13 


8.10 ± 0.04 


UM228, SHOC 11 


SDSS J003218. 60+150014.2 


J0032 




H06 (WHT) 




0.018 


40.30 


7.93 ± 0.03 


SHOC 22 


SDSS J145506. 06+380816.6 


J1455 


H08 


(CAHA - 3.5 


m.) 


0.028 


40.84 


7.94 ± 0.03 


CG 576 


SDSS J150909. 03+454308.8 


J1509 


H08 


(CAHA - 3.5 


m.) 


0.048 


41.26 


8.19 ± 0.03 


CG 642 


SDSS J152817.18+395650.4 


J1528 


H08 


(CAHA - 3.5 


m.) 


0.064 


41.63 


8.17 ± 0.04 




SDSS J154054.31+565138.9 


J1540 


H08 


(CAHA - 3.5 


m.) 


0.011 


39.65 


8.07 ± 0.05 


SHOC 513 


SDSS J161623. 53+470202.3 


J1616 


H08 


(CAHA - 3.5 


m.) 


0.002 


38.93 


8.01 ± 0.03 




SDSS J162410. 11-002202.5 


J1624 




H06 (WHT) 




0.031 


41.48 


8.05 ± 0.02 


SHOC 536 


SDSS J165712. 75+321141.4 


J1657 


H08 


(CAHA - 3.5 


m.) 


0.038 


40.74 


7.99 ± 0.04 




SDSS J172906. 56+565319.4 


J 1729 


H08 


(CAHA - 3.5 


m.) 


0.016 


40.57 


8.08 ± 0.04 


SHOC 575 



SED and the analysis of the Wolf-Rayet features. We also 
describe the photoionization models of the gas, once the 
underlying stellar population has been removed. In section 
4, we discuss the star formation history and the physical 
properties obtained from the models and we compare them 
with observations. In section 5, we summarize our results 
and present our conclusions. Finally, an appendix has been 
added to compare the total abundances and ionization cor- 
rection factors obtained from models with those obtained in 
H06 and H08. 



2 THE SAMPLE OF STUDIED OBJECTS 

The modeled galaxies were selected from the SDSS using the 
implementation of the database in the INAOE Virtual Ob- 
servatory superserveiQ (see H06 and H08 for more details) . 
We selected the brightest nearby narrow emission line galax- 
ies with very strong lines and large equivalent widths of the 
Hq line from the whole SDSS data release available when we 
planned each observing run. Once AGN-like objects were re- 
moved by BPT diagnostic diagrams (Baldwin, Phillips and 
Terlevich, 1981), the sample was restricted to the largest Ha 
flux and highest signal-to-noise ratio objects (Lopez, 2005). 
The HII galaxies studied in this work comprise 10 objects 
belonging to the final list and observed in two independent 
observing runs with different double-arm spectrographs in 
order to improve the spectral range and the signal-to-noise. 
Three out of them were observed with the ISIS spectrograph 
mounted on the William Herschel Telescope (H06) and the 
other 7 were observed with the TWIN instrument in the 3.5 
m. CAHA telescope (H08). In Table [1] we list the galaxies, 
with the names adopted in this work, the spectroscopic red- 
shift as listed in the SDSS catalogue and their extinction- 
corrected Ha luminosities in erg-s^^. Distances have been 
derived from the measured redshifts for cosmological pa- 
rameters: Ho = 73 km S~^ MpC~^, Qmatter = 0.27, ^vaccum 

= 0.73. We also list the oxygen abundance as calculated in 
H06 and H08 and, finally, other names by which the galaxies 
are also known. 

Both observing runs had similar configurations. In 
WHT observations the spectral range goes from 3200 A up 
to 10550 A with a spectral dispersion of 2.5 A/px FWHM 
in the blue arm and 4.8 in the red. In CAHA observations, 
the spectral range covers from 3400 up to 10400 A (with 

^ http://ov.inaoep.mx/ 



a gap between 5700 and 5800 A) with spectral resolutions 
of 3.2 and 7.0 A/px FWHM in the blue and the red part 
respectively. The observations were taken at parallactic an- 
gle. These configurations allow to detect simultaneously the 
bright emission lines from [On] at 3727 A up to [Slll] at 
9532 A, allowing the derivation of five different electron tem- 
peratures (t([On]), t([Oin]), t([Sn]), t([Sin]) and t([Nn]) ) 
in three objects of the sample and four temperatures (the 
same listed above with the exception of t([Nn]) in the other 
seven). The derivation of all these electron temperatures al- 
lows a more precise determination of the thermal structure 
of the gas and, hence, of the ionic abundances of the de- 
tected species in the observed spectra. Our sample galaxies 
have oxygen abundances in the range 7.93 < 12+log(0/H) 
< 8.19 and equivalent widths of H/3 in the range 83 A < 
EW(H/3) < 171 A. 



3 RESULTS 

3.1 Model fitting of the stellar population 

Underlying stellar populations in starburst galaxies have 
several effects in their spectra that can affect the measure- 
ment of emission lines. Firstly, since this population has a 
weight in the galaxy continuum, most of the times the mea- 
surement of the equivalent widths of the emission lines is 
not a trustable estimator of the age of the ionizing stellar 
population (Terlevich et al., 2004). Secondly, the presence 
of a conspicuous underlying stellar population depresses the 
Balmer and Paschen emission lines and do not allow to mea- 
sure their fiuxes with an acceptable accuracy (Di'az, 1988). 
This can alter those properties dependent on the fluxes of 
these recombination lines, like the reddening or the ionic 
abundances. In the case of the helium recombination lines, 
this effect is present too, affecting the derivation of the he- 
lium abundance (Olive & Skillman, 2004). 

We have subtracted from the observed spectra the spec- 
tral energy distribution of the underlying stellar popula- 
tion found by the spectral synthesis code STARLIGHT 
(Cid Fernandes et al. 2004, 2005; Mateus et al. 2005). 
STARLIGHT fits an observed continuum spectral energy 
distribution using a combination of multiple simple stellar 

^ The STARLIGHT project is supported by the Brazilian agen- 
cies CNPq, CAPES and FAPESP and by the France-Brazil 
CAPES/Cofecub program 
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Figure 1. In the solid brown line, optical spectrum of the galaxy J0021 and comparison with the spectral fitting carried out with 
the code STARLIGHT, using spectral libraries from Starburst99, represented by the solid red line. In the lower panel, we show the 
residuals of the fit, which represents the nebular spectrum. 




3700 
Wavelength ( A) 



Figure 2. Upper panel, in black, optical spectrum of the galaxy 
J1624 in the spectral range 3500 - 3912 A, around the Balmer jump 
and the high order Balmer series, together with the fitting (red 
line) performed using the STARLIGHT spectral synthesis code. 
Lower panel, subtraction of this fitting from the WHT spectrum 
in the same spectral range; the solid red line shows the zero level. 



population synthetic spectra (SSPs; also known as instanta- 
neous burst) using a minimization procedure. In order to 
keep the consistency with the photoionization models used 
to model the gas emission, we have chosen for our anal- 
ysis the SSP spectra from the Starburst99 libraries (Lei- 
therer et al., 1999; Vazquez & Leitherer, 2005), based on 
stellar model atmospheres from Smith et al. (2002), Geneva 
evolutionary tracks with high stellar mass loss (Meynet et 
al., 2004), a Kroupa Initial Mass Function (IMF; Kroupa, 
2002) in two intervals (0.1-0.5 and 0.5-100 Mq) with dif- 
ferent exponents (1.3 and 2.3 respectively), the theoretical 
wind model (Leitherer et al., 1992) and a supernova cut-ofT 
of 8 M0. We have fixed the metallicity of the stellar pop- 
ulations to Z = 0.004 (= 1/5 Zq) and Z = 0.008 (= 2/5 
Zq) depending on the closest total oxygen abundance as 
measured using the direct method for each galaxy ( see HOG 
and H08). The STARLIGHT code solves simultaneously the 
ages and relative contributions of the different SSPs and the 
average reddening. The reddening law from Cardelli, Clay- 
ton & Matthews (1989) has been used. Prior to the fitting 
procedure, the spectra were shifted to the rest-frame, and 
re-sampled to a wavelength interval of 1 A in the entire wave- 
length range between 3500 A and 9000 A by interpolation 
conserving flux, as required by the program. Bad pixels and 
emission lines were excluded from the final fits. 

To illustrate the fitting results, we show in the upper 
panel of Figure [T] the spectrum of J0021 with the best fit 
found by STARLIGHT overimposed and, in the lower panel 
of the same figure, the emission line spectrum once the un- 
derlying stellar continuum has been subtracted. The mea- 
sured intensities of the most prominent Balmer emission 
lines in this subtracted spectrum are, in all cases, equal 
within the errors to those reported in H06 and H08, where a 
pseudo-continuum was adopted. It was also shown in those 
works that the measurement of the intensities of the Balmer 
lines is consistent with a multi-gaussian fitting of the lines. 
This method fits independently the emission line and the 
absorption feature, which is visible through their wings in 
some bright Balmer lines. In the upper panel of Figure 2 
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Figure 3. Histogram of the distribution of the percentage of visual Hght from each population in the best solution found by 
STARLIGHT in the spectral fitting of the optical spectrum for each galaxy of the sample. 



we plot in detail the WHT spectrum of J 1624 in the spec- 
tral range 3500-3912 A, around the Balmer jump and the 
high order Balmer series, and in solid red line we show the 
STARLIGHT fitting. In the lower panel of the same figure 
we show the subtraction of the fit from the observed spec- 
trum. We can appreciate that small errors in the fit of the 
underlying stellar population could become a great unquan- 
tifiable error in the emission line fiuxes. Nevertheless, for the 
strongest emission lines, the differences between the mea- 
surements done after the subtraction of the STARLIGHT 
fit and those using the pseudo-continuum approximation are 
well below the observational errors, giving almost the same 
result. It must be highlighted that in this type of objects the 
absorption lines, used to fit the underlying stellar popula- 
tions, are mostly affected by the presence of emission lines. 
There are a few that are not, such as some calcium and 
magnesium absorption lines. To make the STARLIGHT fit 
it is necessary to mask the contaminated lines, thus the fit is 
based on the continuum shape and takes into account only 
a few absorption lines. Hence, it is not surprising that the 
STARLIGHT results are not so good for this kind of ob- 
jects. It must be noted too that the STARLIGHT method- 
ology is very advantageous for statistical and comparative 
studies. Those studies only consider the results derived from 



the strongest emission lines which have small proportional 
errors. However, we must be careful when the aim of our 
work is the detailed and precise study of a reduced sample 
of objects because, as we have shown, we can not quantify 
the error introduced by the fitting. Moreover, these errors 
are higher for the weakest lines, including the auroral ones. 
Regarding helium absorption features the multigaussian ap- 
proximation is not possible in any line, because their wings 
are not visible. Nevertheless, we have checked that the con- 
tribution of the depressed components of the emission lines 
predicted by the STARLIGHT fit are not larger than taking 
into account a slightly lower position of the continuum, as 
it is explained in both HOG and H08. 

/,From this analysis we only take as a valuable input 
information for our tailor-made photoionization models, the 
ratios predicted by STARLIGHT for the ionizing and the un- 
derlying stellar populations in order to correct the equivalent 
width of H/3 and to take it as an estimate of the properties 
of the ionizing stellar population. 

Nevertheless, we can do a comparative analysis of the 
star formation histories of every galaxy as predicted by 
the population synthesis fitting technique. In Figure O we 
show the histogram of the age distribution of the contri- 
bution to the visual hght of each stellar population used 
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Table 2. Properties of the stellar populations as obtained with STARLIGHT for each studied object. In the last two columns we list 
the measured EW(H/3) and its value corrected for contribution to the continuum from the young stellar component only. 



Object ID 


A(V) 


log M, 


log Mion, 


% Mio„. 


EW{H/3) obs. 


EW(H/3) cor 




(mag) 


(Me) 


(Me) 




(A) 


(A) 


J0021 


0.01 


9.40 


6.89 


0.31 


97 


110 


J0032 


0.17 


8.15 


5.53 


0.24 


90 


119 


J1455 


0.09 


8.42 


6.04 


0.42 


133 


146 


J1509 


0.58 


8.08 


6.96 


7.57 


123 


135 


J1528 


0.25 


8.92 


6.81 


0.78 


171 


191 


J1540 


0.42 


7.51 


5.22 


0.52 


122 


124 


J1616 


-0.11 


6.55 


4.50 


0.89 


83 


83 


J1624 


0.60 


8.88 


5.73 


0.07 


101 


107 


J1657 


0.16 


8.45 


6.20 


0.56 


118 


132 


J1729 


0.27 


8.65 


6.07 


0.26 


126 


126 



by STARLIGHT to fit the optical SED for the ten studied 
galaxies As we can see, tlie blue light in these objects seems 
to be largely dominated by a very young stellar population 
(younger than 10 Myr) which is responsible for the ioniza- 
tion of the gas. At same time, there are different underlying 
stellar populations, with a component older than 1 Gyr in 
all cases. This "old" component dominates the total stellar 
mass in all the objects, providing more than 99 % of it, with 
the exception of J1509, for which this percentage is about 92 
%. At any rate, we find very similar age distributions for the 
stellar populations in all the objects pointing to very sim- 
ilar star formation histories with the presence of recursive 
starburst episodes resembling the one observed at present. 

In table [21 we summarize the properties of the best fit- 
ting to the optical spectrum of each galaxy as found by 
STARLIGHT. We list in columns 2 to 5 the visual extinc- 
tion, A(V), the total stellar mass, M*, the mass of the stellar 
population younger than 10 Myr, Mion* , and the percentage 
contribution of this to the total stellar mass, % Mion*- We 
also list in columns 6 and 7 the measured equivalent width 
of H/3, EW(H/?), and the corrected value, once the contribu- 
tion to the continuum by the underlying stellar population 
has been subtracted. As we can see, our procedure yields a 
negative value for the extinction in J1616. Nevertheless, this 
is still consistent with the very low extinction found using 
the Balmer decrements. 

We can check the accuracy of these fittings by compar- 
ing the stellar mass of the ionizing population as derived 
from STARLIGHT, and the values of the same masses as 
derived from the extinction-corrected Ha fluxes, using the 
expression from Diaz (1998): 

log Mio„, = log L{Ha) - 0.86 ■ log EW{HI3) - 32.61 (1) 

This comparison is shown in Figure 2] It can be seen 
that the agreement between both derived values is very good 
for most of the studied objects. Only J0021 shows a higher 
young stellar mass when derived from its Ha flux. Horizon- 
tal error bars in the plot correspond to the underestimate of 
the Ha flux due to the UV dust absorption. We have quan- 
tifled this effect taking the absorption factors derived from 
the photoionization models described in section 3.3 below. 
Nevertheless, taking into account these corrected values of 
L(Ha) in the derivation of M^on* from (1) does not increase 
the ionizing stellar masses by more than 0.3 dex, except in 
the case of J1455, with 0.42 dex. 




4 4,5 5 5,5 6 6,5 7 7,5 8 8,5 



log (M.^^ ,) (from Ha flux) 

Figure 4. Comparison between the ionizing stellar mass as de- 
rived from the extinction-corrected Ha fluxes and from SED fit- 
ting for an age younger than 10 million years. The red solid line 
represents the 1:1 relation. Error bars to the right take into ac- 
count the dust absorption factors obtained from the photoioniza- 
tion models described below. 

3.2 Detection and measurement of Wolf-Rayet 
bumps 

Wolf-Rayet (WR) stars appear in the first stages after the 
main sequence phase of the evolution of massive stars, so 
they are already visible soon after the beggining of a burst 
of star formation (approximately 2 Myr) and for a relatively 
short period of time (about 3 Myr, on average). The exis- 
tence of this WR stellar population can be very useful for 
the study and characterization of the ionizing stellar pop- 
ulation in starburst galaxies (e.g. Perez-Montero & Diaz, 
2007). The strength of the stellar winds produced by these 
stars can be sometimes measured in the integrated spec- 
tra of the starburst galaxies, which are then identified as 
WR galaxies (Conti, 1991). The presence of the emission 
broad features ("bumps") produced by WR stars are there- 
fore associated to the existence of processes of intense star 
formation. These WR features are the blue bump, centered 
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Table 3. Observed extinction corrected H/3 intensities and properties of the WR bumps for the objects of the sample. We list equivalent 
widths and relative intensities, normalized to 100 times the H/3 value, measured at 4650 A (blue bump, bb) and at 5808 A (red bump, rb). 
The corrected values take into account the contribution of the underlying stellar population and the nebular emission to the continuum 
and the dust absorption of the UV predicted by the photoionization models in each object. 



Object ID 


I(H/3) 


-EW(bb) 


-EW(bb),o, 






-EW(rb) 


-EW(rb),<,, 


I(rb)" 


I(i-b)S„, 




(erg ■ s cm ^) 


(A) 


(A) 






(A) 


(A) 






J0021 


2.47 ■ 10"^"* 


3.9±0.9 


5.8 


1.4±0.3 


0.8 










J0032 


1.23 ■ 10"^'' 


S.litl.O 


6.7 


4.9±1.0 


3.2 


1.3±0.6 


1.9 


0.9±0.4 


0.6 


J1455 


1.49 ■ 10-14 


















J1509 


1.35 ■ 10"" 


9.1±2.0 


12.0 


8.3±1.8 


5.3 


1.6±0.5 


2.2 


1.1±0.4 


0.7 


J1528 


1.73 ■ 10"" 


















J1540 


5.30 ■ 10-15 


10.3±2.8 


16.5 


9.8±2.6 


9.5 










J1616 


1.36 ■ 10-1'' 


5.5±3.8 


6.6 


6.9±4.8 


3.3 










J1624 


5.03 ■ 10-14 


7.4±1.5 


13.1 


2.6±0.6 


1.3 


2.8±0.9 


3.8 


0.8±0.3 


0.4 


J1657 


6.30 ■ 10-15 


















J 1729 


2.40 ■ 10-14 


6.5±2.1 


9.4 


6.5±2.1 


3.4 











in units of 100-I(H/3) 



at 4650 A and produced mainly by broad emission lines of 
Nv at 4605, 4620 A, Niii 4634, 4640 A, Ciii/iv 4650, 4658 
A and Hell at 4686 A, and the red bump, usually fainter, 
centered at 5808 A and emitted mainly by Cm. 

We have detected the blue bump in seven of the ten 
galaxies of our sample, while the red bump has been de- 
tected in only three of these seven galaxies. All the three 
galaxies observed with WHT had already been identified as 
WR galaxies in H06 and have also been included in the cata- 
logue of WR galaxies of the SDSS (Brinchmann et al. , 2008) . 
The other four galaxies with a detection of WR features, ob- 
served with the 3.5 CAHA telescope are not listed in that 
catalogue. An example of the WR emission can be seen in 
Figure 6 of H06 for the J 1624 galaxy. 

We have measured the bumps in the following way. First 
we have adopted a continuum shape and we have measured 
the broad emission over this continuum at the wavelengths 
of the blue bump, at 4650 A, and the red bump, at 5808 
A. Finally, we have subtracted the emission of the narrower 
Unes not emitted by the WR stars winds, which are [Fein] 
4658 A, the narrow component of Hen 4686 A, Hel + [Ariv] 
4711 A and [Ariv] 4740 A in the blue bump. The equiv- 
alent widths and the deredenned relative intensities of the 
blue and red bumps are shown in Table [3] Both have been 
corrected in order to study the properties of the ionizing 
stellar population. In the case of EWs, we have subtracted 
the contribution of the underlying stellar population as de- 
rived from STARLIGHT and from the nebular continuum 
as derived from the photoionization models described below. 
In the case of the relative intensities, we have taken into ac- 
count the dust absorbed component of H/3 derived from the 
same photoionization models for each object. 

3.3 Photoionization models of the ionized gas 

Photoionization models have been calculated using the code 
CLOUDY V. 06.02c (Ferland et al., 1998) with the same 
atomic coefficients used to derive ionic abundances in H06 
and H08 and the recent dielectronic recombination rate co- 
efficients from Badnell (2006). We have taken the same Star- 
burst99 stellar libraries used by STARLIGHT to fit the ob- 
served SEDs. We have assumed different initial input param- 



eters and then, we have followed an iterative method to fit 
the relative intensities of the strongest emission lines in the 
integrated spectra, along with other observable properties. 

We have assumed a constant star formation history, 
with a star formation rate calculated according to the mea- 
sured and corrected number of ionizing photons. The age of 
this burst of star formation has been limited to the range 1- 
10 Myr in order to reproduce the corrected equivalent width 
of H/3 once the older stellar populations have been removed, 
except in the case of J1540 for which a burst duration of 
15.85 Myr has been considered in order to reproduce the ob- 
served EW(H/3). Total abundances have been set to match 
the values derived in H06 and H08, including He, O, N, 
S, Ne, Ar and Fe. The abundances of the rest of the ele- 
ments have been scaled to the measured oxygen abundance 
of each galaxy and taking as reference the photosphere so- 
lar values given by Asplund et al. (2005). The inner radii of 
the ionized regions have been chosen to reproduce a thick 
shell in all cases. This kind of geometry reproduces best 
the ionization structure of oxygen and sulphur simultane- 
ously (Perez-Montero & Di'az, 2007). We have assumed a 
constant density equal to that derived from the measured 
[Sll] emission line ratio (~ 100 particles per cm'' in all the 
objects). We have kept as free parameters the filling factor 
and the amount of dust. The inclusion of dust allows to fit 
correctly the measured electron temperatures because it is 
a fundamental component to reproduce the thermal balance 
in the gas. The heating of the dust can affect the equilibrium 
between the cooling and heating of the gas and, thus, the 
electron temperature inside the nebula. We have assumed 
the default grain properties given by Cloudy 06, which has 
essentially the properties of the interstellar medium and fol- 
lows a MRN (Mathis, Rumpl & Norsieck 1997) grains size 
distribution. 

After some first estimative guesses of the UV absorp- 
tion factor due to the internal extinction, fd, we have re- 
calculated the age of the cluster and the number of ionizing 
photons to match the observed values. Finally, we have ad- 
justed the chemical abundances to fit the relative intensities 
of the emission lines. 

The average number of models needed to fit the observ- 
able properties of the galaxies with a reasonable accuracy 
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is about 50 for each galaxy. In Table [4] we show the inten- 
sities of the most representative emission lines as measured 
in HOG and H08 compared to the values obtained with the 
best model for each object. The agreement in most cases is 
excellent, even for the faint auroral lines. 

In Table O we show the different line temperatures ob- 
tained from the photoionization models compared to those 
derived in H06 and H08. In the same table, we show the 
elemental ionic abundances derived using the correspond- 
ing electron temperature representative of the region where 
the ion is formed: t([Olll]) for 0^+, ^e^+ and Fe^+; t([Slll]) 
for and Ar^+i t([Oii]) for 0+ and t([Nn]) for N+ when 
t([Nn]) is available or t([On]) in the other cases. These abun- 
dances are compared with the ionic abundances obtained 
from the corresponding models. We also show the tempera- 
ture from the oxygen recombination lines predicted by our 
models in the high excitation region and the temperature 
fluctuations ijp' , Peimbert, 1967) obtained from the differ- 
ence between the temperatures derived from recombination 
and coUisional emission lines. It can be seen that the pre- 
dicted fluctuations are practically zero and, therefore, no 
ADFs are expected according to our models. The compari- 
son between the model predictions and the temperature fluc- 
tuations inferred from the observed values of t([Olll]) and 
the Balmer jump temperature for the three objects observed 
with the WHT gives compatible results, with the exception 
of J0032, whose is higher as obtained from the observa- 
tions. 



4 DISCUSSION 

4.1 Properties of the ionizing population 

In Table|B]we list the main observed properties of the young 
stellar population responsible for the ionization of the sur- 
rounding gas and their comparison with the same values 
as predicted by the corresponding tailor-made models. Col- 
umn 2 lists the age of the burst of continuous star formation 
adopted for each model in Myr. Column 3 gives the metal- 
licity {Z) as derived from the relative oxygen abundance 
measured in the gas-phase and column 4 the adopted Z of 
the stellar ionizing cluster, which is in each case the value 
closest to the oxygen abundance. Although the hypothesis of 
equal metallicity for the gas and the ionizing stars is not well 
justified and the available metallicities in the stellar libraries 
are still poor to some extent, the differences in metallicity 
among models are not crucial to improve their accuracy. 
Column 5 in the table lists the number of ionizing photons, 
Qobs(H), as derived from the observed luminosity of Ha us- 
ing the following expression (see, for example, Osterbrock, 
1989): 



,(//) = 7.35 • 10" -LiHa) 



(2) 



with L(Ha) in erg ■ s~^. 

This Qabs (H) should be compared to the the number of 
photons produced by the model cluster once divided by the 
dust absorption factor, fd, predicted by the models, which 
is listed in column 6. The total number of ionizing photons 
produced by the ionizing cluster is then: 



where both Q(H) are expressed in photons ■ s~^. Therefore, 
in order to be compared with the observations we list in col- 
umn 7 the ratio Q{H)/ fd predicted by the models. Finally, 
columns 7 and 8 give the equivalent width of H/3, EW(H/3), 
measured on the spectra once the underlying stellar popu- 
lations have been removed. This is to be compared with the 
corresponding EW(H/3) value given by the models, which al- 
ready takes into account the dust absorption factor and the 
nebular continuum emission. In the lower panel of the table 
we show in columns 2 and 3 the logarithm of the ionization 
parameter, log U, calculated using the following expression 
derived by Diaz et al. (1991) from the emission lines of [Sll] 
and [Sin]: 



\ogU 



MogJ^^^MMZ5i4_2.99 



(4) 



Q{H) = fd-Qobs{H) 



(3) 



/([5'm]9069,9532A) 

which ofi^ers a reliable estimation of the ionization parameter 
in nebulae with a non-plane parallel geometry. This value is 
compared with the corresponding log U obtained from our 
models. The logarithm of the filling factor (log e) and the 
adopted inner radius (Ro) of the ionized region are listed in 
columns 4 and 5 while columns 6 and 7 give the external 
radius (Rs) from the models as compared with the radius 
derived from the spatial profile of Ha in the long slit. Finally, 
the last three columns give the dust-to-gas ratio in each 
model and the visual extinction, A(V), as derived from the 
observed Balmer decrements and compared with the values 
obtained from the models. 

As in the case of the emission lines, the agreement be- 
tween the derived properties and those obtained from the 
models is good with the exception of the outer radii which 
are much larger in the observations by factors ranging from 
2.2 (J1616) to 15.5 (J1540). Possible explanations for these 
discrepancies could be either a very low density layer of dif- 
fuse gas in the outer regions of the ionized gas or lower filling 
factors as compared with the values found in our models. We 
find as well some discrepancies between the extinction values 
found in our models and those derived using Balmer decre- 
ments. This could be symptomatic of very complex struc- 
tures of dust in the ionized gas region, not easy to reproduce 
in the models. 

Regarding the age of the burst, the equivalent width of 
an intense Balmer emission line, like H/3, is generally associ- 
ated to the age and properties of the ionizing population. It 
has long been accepted that the low values of EW(H/3) found 
in starburst galaxies as compared with the predictions from 
Population Synthesis Models are due to the presence of an 
underlying stellar population that enhances the continuum 
light without significantly contributing to the ionization of 
the gas. Nevertheless, as can be seen in Table [2]the subtrac- 
tion of this population from the integrated spectra does not 
imply significant corrections to EW(H/?) (between 1.5 % in 
the case of J1540 and J1729 and 24 % in the case of J0032). 

However, the excitation conditions found in these ob- 
jects, are only reachable with the energy supplied by very 
young and massive O stars, whose associated EW(H/3) is 
four or five times larger than the observed values. This dif- 
ference can be partially explained by a certain dust absorp- 
tion of UV photons. The dust absorption factors obtained 
in our models range between 1.5 and 2, as we can see in 
Table |6l These help to put the modeled EW(H/3) values in 
agreement with the observed ones. 
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Table 4. Comparison between the observed intensities of the most representative emission Hncs in terms of 100 ■ F(H/3) and the 
reproduced intensities by the best model obtained for each object. 



A 


(A) 


J0021 




J0032 




,71455 




J1509 




J1528 








Observed 


Model 


Observed 


Model 


Observed 


Model 


Observed 


Model 


Observed 


Model 


3727 


[On] 


163.4 ± 1.9 


162.9 


157.3 ± 1.4 


157.9 


111.5 ± 1.6 


117.5 


153.2 ± 1.8 


163.6 


228.8 ± 2.9 


230.5 


3868 


[Nclll] 


38.8 ± 0.8 


40.6 


37.2 ± 1.0 


39.2 


47.9 ± 1.1 


50.1 


35.0 ± 1.1 


34.1 


47.2 ± 1.8 


46.6 


4363 


[Oiii] 


5.6 ± 0.3 


5.7 


6.2 ± 0.3 


6.0 


10.2 ± 0.4 


10.3 


4.2 ± 0.2 


4.1 


5.0 ± 0.2 


5.0 


4658 


[Felll] 


1.0 ± 0.1 


1.1 


0.9 ± 0.1 


0.9 


- 


1.0 


1.1 ± 0.1 


0.9 


1.3 ± 0.2 


1.4 


4959 


[Oiii] 


153.2 ± 1.3 


140.9 


158.2 ± 1.3 


156.1 


204.6 ± 1.3 


192.5 


167.5 ± 1.5 


154.0 


165.6 ± 1.7 


162.2 


5007 


[Oiii] 


433.4 ± 2.6 


424.0 


460.7 ± 2.4 


469.9 


613.6 ± 3.4 


579.4 


499.4 ± 1.5 


403.4 


489.3 ± 2.92 


488.1 


5876 


Hei 


12.7 ± 0.7 


IS.l 


12.4 ± 0.6 


11.4 


11.4 ± 0.35 


10.2 


12.68 ± 0.59 


12.6 


12.27 ± 0.36 


11.8 


6312 


[Siii] 


1.2 ± 0.1 


1.6 


2.3 ± 0.1 


2.1 


1.6 ± 0.1 


1.8 


1.6 ± 0.1 


2.0 


1.7 ± 0.1 


1.6 


6548 


[N II] 


— 


9.5 


3.7 ± 0.2 


3.6 


2.7 ± 0.2 


3.0 


5.4 ± 0.3 


5.0 


7.2 ± 0.5 


7.0 


6584 


[Nil] 


26.0 ± 0.9 


28.1 


11.2 ± 0.5 


11.2 


7.9 ± 0.2 


8.8 


13.9 ± 0.4 


14.6 


19.6 ± 0.4 


20.6 


6717 


[Sll ] 


13.6 ± 0.5 


14.5 


17.3 ± 0.5 


10.9 


10.0 ± 0.3 


9.6 


19.7 ± 2.1 


18.3 


19.2 ± 0.8 


18. 4 


6731 


[Sii] 


10.7 ± 0.5 


11.0 


12.5 ± 0.3 


12.7 


7.9 ± 0.2 


7.2 


14.9 ± 0.4 


13.9 


14.2 ± 0.6 


14.0 


7136 


[Arm] 


6.2 ± 0.2 


6.5 


9.2 ± 0.3 


8.2 


6.62 ± 0.3 


7.2 


9.82 ± 0.29 


9.8 





8.9 


7319 


[On] 


2.0 ± 0.1 


2.4 


2.5 ± 0.1 


2.4 


1.8 ± 0.1 


1.9 


2.0 ± 0.1 


2.2 


3.0 ± 0.2 


3.1 


7330 


[On] 


1.3 ± 0.1 


2.0 


1.9 ± 0.1 


1.9 


1.4 ± 0.1 


1.5 


1.7 ± 0.1 


1.8 


2.4 ± 0.1 


2.6 


9069 


[Sin] 


11.9 ± 0.6 


15.0 


21.7 ± 1.0 


19.9 


11.5 ± 0.6 


14.1 


25.5 ± 1.2 


23.9 


16.9 ± 1.3 


17.3 


9532 


[Sin] 


23.9 ± 0.9 


37. 2 


44.5 ± 3.8 


47.8 


32.8 ± 1.5 


34.9 


51.8 ± 2.6 


59.3 


40.5 ± 2.9 


43.0 


A (A) 


.J1540 




.11616 




.J1624 




J1657 




J1729 








Observed 


Model 


Observed 


Model 


Observed 


Model 


Observed 


Model 


Observed 


Model 


3727 


[On] 


217.9 ± 2.6 


215.3 


84.9 ± 2.0 


84. 1 


147.1 ± 2.2 


157.6 


188.3 ± 2.3 


203.2 


176.2 ± 2.4 


I64.4 


3868 


[Nclll] 


21.4 ± 0.8 


21.8 


41.1 ± 1.7 


44.0 


42.8 ± 1.1 


41.4 


32.6 ± 1.3 


33.6 


47.9 ± 1.7 


50.0 


4363 


[Oiii] 


2.9 ± 0.2 


2.9 


8.5 ± 0.3 


9.1 


7.0 ± 0.2 


7.1 


5.2 ± 0.2 


5.2 


6.6 ± 0.3 


6.5 


4658 


[Fein] 


0.8 ± 0.1 


0.8 




6.9 


0.9 ± 0.1 


0.7 


1.1 ± 0.2 


1.1 


0.9 ± 0.1 


0.9 


4959 


[Oiii] 


104.8 ± 0.8 


115.3 


204.9 ± 1.5 


204.8 


193.8 ± 1.4 


177.5 


143.3 ± 1.3 


130.0 


171.0 ± 1.5 


184. 3 


5007 


[Oiii] 


309.4 ± 1.9 


347.2 


615.2 ± 3.7 


616.5 


564.2 ± 2.8 


534.4 


430.8 ± 2.4 


391.4 


515.4 ± 4.2 


554.7 


5876 


Hoi 


11.6 ± 0.4 


12.8 


10.6 ± 0.8 


11.2 


13.5 ± 0.3 


14.0 


11.2 ± 0.4 


11.4 


12.6 ± 0.4 


13.1 


6312 


[Sin] 


1.4 ± 0.1 


1.8 


1.9 ± 0.1 


1.9 


1.8 ± 0.1 


1.7 


2.0 ± 0.1 


1.8 


1.8 ± 0.1 


1.8 


6548 


[Nil] 


7.4 ± 0.3 


5.5 


1.6 ± 0.1 


1.3 


3.4 ± 0.2 


3.0 


4.6 ± 0.2 


4.9 


7.7 ± 0.2 


7.6 


6584 


[Nil] 


21.2 ± 0.7 


16. 4 


4.3 ± 0.2 


4.0 


9.3 ± 0.2 


8.9 


14.3 ± 0.5 


14.6 


21.9 ± 0.5 


22.6 


6717 


[Sii] 


26.1 ± 0.5 


25.2 


7.7 ± 0.2 


7.9 


13.6 ± 0.3 


12.6 


22.1 ± 0.6 


19.6 


12.9 ± 0.3 


14.4 


6731 


[Sii] 


19.2 ± 0.5 


19.1 


5.8 ± 0.2 


6.0 


9.9 ± 0.3 


9.5 


16.0 ± 0.4 


14.8 


10.1 ± 0.3 


10.9 


7136 


[Ann] 


8.9 ± 0.5 


8.8 


7.4 ± 0.4 


7.1 


8.0 ± 0.3 


8.0 


7.2 ± 0.3 


6.3 


8.6 ± 0.4 


9.5 


7319 


[On] 


2.7 ± 0.2 


2.9 


1.4 ± 0.1 


1.3 


2.1 ± 0.1 


2.4 


3.0 ± 0.2 


3.0 


2.3 ± 0.1 


2.4 


7330 


[On] 


2.2 ± 0.1 


2.3 


0.9 ± 0.1 


1.1 


1.8 ± 0.1 


1.9 


2.1 ± 0.1 


2.4 


1.9 ± 0.1 


1.9 


9069 


[Sin] 


21.0 ± 0.6 


22.0 


16.5 ± 0.6 


15.9 


7.0 ± 0.5 


15.4 


14.0 ± 1.0 


16.6 


20.9 ± 1.0 


18.1 


9532 


[Sin] 


53.3 ± 3.3 


54.7 


40.1 ± 1.5 


39.5 


40.1 ± 1.9 


38.2 


36.7 ± 2.6 


41.2 


47.2 ± 2.5 


45.0 
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Figure 5. Left panels: relation between the equivalent width of H/3 and the age of the corresponding synthesis population model 
for metallicities of 0.004 and 0.008 as predicted by Starburst99. Crosses represent the observed values measured on the young 
stellar continua and multiplied by the corresponding absorption factor obtained in each tailor-made model. The assigned age to 
each object is that predicted by the same models. Therefore, the J1540 point does not appear as the fitted age is older than 10 Myr. 
Right panels: relation between the ratio of helium to hydrogen ionizing photons and the age of the ionizing cluster. In all panels, 
the blacli solid line represents an instantaneous star formation and the red dashed line, a continuous star formation as predicted 
by Starburst99. As in the previous case, crosses represent the values obtained from the models for each object, with the exception 
of J1540. 



In the left panels of Figure [5] we show the evolution of 
EW(H/3) with time for the two metallicities closest to those 
in our sample for both instantaneous (black solid line) and 
continuous (dashed red line) star formation laws as obtained 
from Starburst99 models. In the right panels of the same fig- 
ure we show the evolution of the ratio of helium to hydrogen 



ionizing photons, Q(He)/Q(H) according to the same Star- 
burst99 models. This ratio, that can be considered a good 
estimate of the ionizing power of the cluster, falls dramati- 
cally for an instantaneous star formation law once the WR 
phase has finished (~ 5 Myr), while it maintains a constant 
high value in the case of a continuous star formation. We 
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Table 5. Comparison between the derived electron temperatures and ionic abundances derived in H06 and H08, and the values obtained 
from the best tailor-made photoionization model for each object. 





.T0021 


Ad odd 


,10032 


Model 


,11455 


Model 


,11509 


Model 


J1528 


Model 


T([Olll) (K) 
12+log(0+/H+) 


10300 ± 200 
7.73 ± 0.06 


1230S 
l.U 


13500 ± 400 
7.15 ± 0.06 


12388 
7.43 


13300 ± 700 
7.16 ± 0.09 


13062 
7.20 


11800 ± 500 
7.48 ± 0.08 


11561 
7.57 


11700 ± 500 
7.67 ± 0.09 


11695 
7.68 


Ti' rOrrrl ^ f 

1 (_[U111J_) {IS-) 

T,([Ollll) (K) 

12+log(02+/H+) 
12+log(No^+ /H+) 
12+log(Fc^+/H+) 


004+""" 
7.86 ± 0.02 
7.27 ± 0.03 
5.50 ± 0.05 


12882 
13000 
0.002 
7.81 
7.16 
5.39 


0.066±0.026 
7.86 ± 0.02 
7.30 ± 0.06 
5.42 ± 0.06 


12531 
12900 
0.001 
7.88 
7.18 
5. S3 


1 4nnn 4- 9nn 

7.87 ± 0.02 
7.20 ± 0.03 
4.81 ± 0.08 


14700 
0. 005 
7.83 
7.12 
5.27 


1 ncinn ^ i nn 

8.10 ± 0.02 
7.44 ± 0.03 
5.71 ± 0.06 


1 1066 
11200 
0. 002 
8.03 
7.28 
5.45 


1 1 finn 4- 1 nn 

8.00 ± 0.02 
7.47 ± 0.04 
5.69 ± 0.07 


11800 
0.001 
7.98 
7.34 
5.58 


T([Nll]) (K) 
12+log(N+/H+) 


11900 ± 500 
6.68 ± 0.04 


12078 
6.50 


6.03 


12162 
6.09 


5.90 ± 0.06 


12647 
5.94 


6.28 ± 0.06 


11482 
6.30 


6.43 ± 0.06 


11561 
6.42 


T([Sll]) (K) 
12+log(S+/H+) 


8600 ± 600 
5.93 ± 0.11 


11776 
5.56 


10300 ± 500 
5.80 ± 0.06 


11885 
5.63 


13100 ± 1100 
5.36 ± 0.06 


12388 
5.32 


8900 ± 700 
6.02 ± 0.10 


11220 
5.74 


9900 ± 700 
5.89 ± 0.10 


11246 
5.71 


12+log(s2+/H+) 
12+log(Ar^+/H+) 


13100 ± 500 
5.91 ± 0.05 
5.50 ± 0.05 


12706 
6.14 
5.52 


13600 ± 700 
6.16 ± 0.06 
5.65 ± 0.06 


12556 
6.26 
5.62 


13700 ± 500 
5.98 ± 0.05 
5.50 ± 0.05 


I4O6O 
6.05 
5.49 


10200 ± 400 
6.44 ± 0.05 
5.94 ± 0.05 


11194 
6.44 
5.81 


12100 ± 600 
6.18 ± 0.07 
5.73 ± 0.09 


11668 
6.26 
5.72 


log(He+/H+) 


-1.05 ± 0.03 


-0.99 


-1.06 ± 0.04 


-1.07 


-1.05 ± 0.05 


-1.11 


-1.00 ± 0.05 


-1.01 


-1.02 ± 0.05 


-1.03 




.J1540 


Ivl ocLgI 


.11616 


Model 


,11624 


Model 


.11657 


Model 


.71729 


Ivl odet 


T([Olll) (K) 
12+log(0+/H+) 


11500 ± 600 
7.67 ± 0.09 


11402 

7.71 


12900 ± 900 
7.07 ± 0.12 


12853 
7.09 


13100 ± 300 
7.18 ± 0.06 


12369 
7.42 


13300 ± 700 
7.37 ± 0.09 


12246 
7.54 


11600 ± 400 
7.57 ± 0.07 


12162 
7.48 


T,([Ollll) (K) 

12+log(02+/H+) 
12+log(Ne^+ /H+) 
12+log(FG2+/H+) 


7.84 ± 0.02 
7.17 ± 0.04 
5.52 ± 0.09 


1 0839 
11000 
0.004 
7. 03 
7.11 
5.40 


1 "^nnn -t- 1 nn 

7.96 ± 0.02 
7.24 ± 0.03 


13428 
13600 
0. 001 

7. as 

7.14 

6.15 


1 9Ann 4- 1 nn 

001+°"^'' 
"■"'"-0.001 

7.98 ± 0.02 

7.33 ± 0.03 

5.48 ± 0.05 


128x53 

13100 

0.002 

7.92 

7.16 

5.28 


1 9'^nn 4- i nn 

7.87 ± 0.02 
7.22 ± 0.04 
5.49 ± 0.08 


12794 

12900 
0.003 
7.78 
7.09 
5.41 


1 9finn 4- 9nn 

7.92 ± 0.02 
7.35 ± 0.04 
5.44 ± 0.08 


12246 
12400 
0.002 
7.99 
7.31 
5.37 


T([Nll]) (A') 
12+log(N+/H+) 


6.47 ± 0.07 


11376 
6.36 


5.67 ± 0.09 


12503 
5.63 


14200 ± 800 
5.94 ± 0.06 


12 134 
6.00 


6.15 ± 0.06 


11995 
6.23 


14000 ± 900 
6.30 ± 0.07 


11967 
6.43 


T([Sll]) (if) 
12+log(S+/H+) 


8500 ± 500 
6.19 ± 0.08 


11143 
5.88 


12100 ± 1200 
5.30 ± 0.10 


12303 
5.26 


10400 ± 700 
5.69 ± 0.08 


11858 
5.49 


8800 ± 500 
6.07 ± 0.08 


11668 
5.70 


8200 ± 600 
5.95 ± 0.10 


1 1 695 
5.58 


T([Slll]) (K) 

12+log(S2+/H+) 

12+log(Ar2+/H+) 


9700 ± 400 
6.47 ± 0.06 
5.95 ± 0.06 


11015 
0.41 
5.78 


12900 ± 700 
6.13 ± 0.05 
5.59 ± 0.06 


IS 366 
6.14 
5. 53 


12600 ± 400 
6.14 ± 0.04 
5.66 ± 0.04 


12706 
6.15 
5.61 


14500 ± 800 
6.00 ± 0.06 
5.49 ± 0.06 


12589 
6.19 
5.51 


11300 ± 500 
6.31 ± 0.06 
5.78 ± 0.06 


12218 
6.25 
5.72 


log(Ho+/H+) 


-1.07 ± 0.03 


-1.07 


-1.08 ± 0.02 


-1.08 


-1.01 ± 0.03 


-1.08 


-1.08 ± 0.06 


-1.08 


-1.05 ± 0.05 


-0.99 



Table 6. Comparison between some properties obtained from the best tailor-made model and derived from the observations for each 
object, including the age of the ionizing cluster, metallicity of the gas and the adopted metallicity of the cluster, dust absorption factor, 
number of ionizing photons, equivalent width of H/3, ionization parameter, filling factor, inner and outer radii, dust-to-gas ratio and visual 
extinction. 



ID Object 


Age (Myr) 


Z 




log Q„6,(H) (s-^) 


Abs. factor log 


! Q(H)//d 


EW(H/3) (A) 






Model 


Observed 


Model 


Observed 


Model 


Model 


Observed 


Model 




J0021 


8.9 


0.0055 ± 0.006 


0.004 


54.00 ± 0.01 


1 . 763 


53.98 


110 ± 2 


112 




J0032 


8.9 


0.0037 ± 0.003 


0.004 


52.17 ± 0.01 


1.508 


52.22 


119 ± 2 


121 




J1455 


4.5 


0.0038 ± 0.004 


0.004 


52.70 ± 0.01 


2.609 


52.58 


146 ± 2 


149 




J 1509 


7.0 


0.0068 ± 0.007 


0.008 


53.13 ± 0.01 


1.500 


53.13 


135 ± 2 


132 




,11528 


5.5 


0.0065 ± 0.007 


0.008 


53.50 ± 0.01 


1.581 


53.48 


191 ± 2 


182 




J 1540 


15.8 


0.0045 ± 0.005 


0.004 


51.51 ± 0.01 


1.027 


51.69 


124 ± 2 


125 




J1616 


10.0 


0.0051 ± 0.004 


0.004 


50.80 ± 0.01 


2.107 


50.78 


83 ± 2 


97 




J1624 


7.9 


0.0049 ± 0.003 


0.004 


53.34 ± 0.01 


1.969 


53.41 


107 ± 2 


111 




J1657 


7.9 


0.0043 ± 0.004 


0.004 


52.60 ± 0.01 


1.621 


52.59 


132 ± 2 


127 




,11729 


7.9 


0.0053 ± 0.005 


0.004 


52.43 ± 0.01 


1.928 


52.30 


126 ± 2 


119 




ID Object 


lof 


5 U 


Fil. fac. log e 


Rq (pc) 


Rs (pc) 




Dust-to-gas 


A(V) 






Observed 


Model 


Model 


Model 


Observed 


Model 




Observed 


Model 


J0021 


-2.71 ± 0.04 


-2.66 


-2.99 


22.7 


5100 ± 200 


1541 


8.76 ■ 10^" 


1.09 ± 0.02 


0.37 


J0032 


-2.41 ± 0.05 


-2.61 


-1.94 


28.5 


1400 ± 200 


177 


5.24 ■ 10^^ 


0.24 ± 0.02 


0.27 


J1455 


-2.33 ± 0.03 


-2.29 


-1.95 


2.9 


2900 ± 600 


240 


9.24 ■ 10~^ 


0.28 ± 0.02 


0.6S 


J1509 


-2.40 ± 0.06 


-2.52 


-2.30 


2.9 


3900 ± 100 


449 


4.19 ■ 10~^ 


0.15 ± 0.02 


0.25 


J1528 


-2.60 ± 0.06 


-2.77 


-2.84 


181.5 


3900 ± 200 


938 


10.58 ■ 10~^ 


0.09 ± 0.02 


O.SO 


J1540 


-2.63 ± 0.04 


-2.84 


-1.62 


67.7 


1200 ± 100 


101 


1.83 . 10"'^ 


0.47 ± 0.02 


0.00 


J1616 


-1.95 ± 0.04 


-2.13 


-0.67 


5.6 


340 ± 70 


22 


5.46 ■ 10^^ 


0.04 ± 0.02 


0.50 


J1624 


-2.48 ± 0.04 


-2.49 


-2.52 


4.8 


2900 ± 200 


684 


7.44 . 10^^^ 


1.12 ± 0.02 


0.42 


J1657 


-2.78 ± 0.05 


-2.84 


-2.50 


15.9 


4000 ± 1000 


374 


10.90 ■ 10^^ 


0.11 ± 0.02 


0.32 


J1729 


-2.20 ± 0.04 


-2.49 


-1.92 


4.4 


1500 ± 100 


182 


6.57 . 10^^ 


0.06 ± 0.02 


0.S7 
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can compare these diagrams with the observed EW(H/3), 
once the underlying stellar population has been removed 
and the dust absorption factor has been taken into account. 
They appear as crosses in Figure [5] and lie in a range which 
corresponds to the asymptotic value of the continuous star 
formation line. In the case of the evolution of Q(He)/Q(H), 
crosses represent the values derived from the tailor-made 
models. 

Although it is true that both EW(H/3) and Q(He) /Q(H) 
could be consistent with an instantaneous burst at present 
in the WR phase, statistically it is not probable that all the 
objects of the sample present such similar properties due to 
their being in the same evolutionary stage. At the same time, 
all the calculated photoionization models in our work require 
a great amount of O and B stars younger than the WR phase 
in order to reach the ionization conditions responsible for the 
observed emission lines. We have therefore assumed for all 
the studied objects a continuous star formation history for 
the ionizing source, with a duration longer in all cases than 
5 Myr. 

4.2 Derived properties of the WR population 

Seven out of the ten observed galaxies have detectable WR 
features in their spectra. Five of them were already included 
in the catalog by Kniazev et al. (2004; SHOC) with a flag in- 
dicating the detection of both the blue and red bumps from 
WR stars. In our spectra, the red bump is observed in only 
three galaxies. The SHOC catalog, however, does not pro- 
vide actual measurements for the features. In a more recent 
paper, Brinchmann et al. (2008) list only the "blue bump" 
luminosities and equivalent widths of three galaxies from the 
present work: J002101 (SHOC 11), J003218 (SHOC 22) and 
J162410 (SHOC 536). The direct comparison between the 
equivalent widths of the blue bump measured by these au- 
thors inside the 3-arcsec fiber spectra of the SDSS and our 
long-slit observations yields discrepancy factors for J0021, 
J0032 and J1624 of 1.7, 2.0 and 1.1, respectively, larger in 
our WHT observations. This could be a consequence of the 
smaller aperture of our long-slit observations which includes 
a lower contribution to the continua from extended ionized 
gas and hence produce larger equivalent widths. 

Table [7] shows the extinction corrected luminosities of 
the measured WR bumps. We show in the same table the 
number of O and WR stars as derived from the extinction 
corrected Ha and WR blue bump luminosities, respectively. 
In the case of O stars, we have taken into account as well the 
dust absorption factor, fd- The total number of stars have 
been estimated by comparison with the properties of the 
Starburst99 models at the age derived from the photoion- 
ization models described above. For those galaxies with a 
red bump detection we have used the same libraries to de- 
rive the ratio between WC and WN stars, which is also listed 
in Table El 

We can compare the relative intensities and equivalent 
widths of the WR bumps detected in our objects with the 
predictions of Starburst99 populations synthesis models as 
a function of the metallicity, age and star formation law of 
the cluster. To make this comparison we have taken the cor- 
rected values listed in Table 3: in the case of EWs from the 
effect of the underlying stellar populations and the contribu- 
tion of the nebular continuum and in the case of the relative 



intensities from the effect of the dust absorption factor to 
the hydrogen recombination lines. In Figure |6] we show the 
predicted equivalent widths and intensities relative to H/3 of 
both the blue and red bumps as a function of the age of the 
cluster for metallicities Z — 0.004 (O.2-Z0), left panels, and 
Z = 0.008 (O.4-Z0), right panels. The black solid lines rep- 
resent the evolution of an instantaneous burst, while the red 
dashed lines do for a continuous star formation history with 
a constant star formation rate. Crosses represent the cor- 
rected values as listed in Table 3 with their corresponding 
error bars. As we can see, the models predict the bright- 
est and most prominent features for the highest metallicity, 
agreeing with the stellar model atmospheres for WR stars 
(Crowther, 2007 and references therein). At the same time, 
we can appreciate that in the instantaneous star formation 
case the WR features appear during a time interval between 
2 and 5 Myr reaching higher intensities than in the case of 
continuous star formation. In this latter scenario, the WR 
features appear at the same age and, despite reaching lower 
intensities, they converge to a non-zero value at older ages. 
For an average metallicity between the two cases of study 
these values are EW ~ 2 A and 0.013-I(H/J) for the blue 
bump and EW ~ 1 A and 0.004-I(H^) for the red one. By 
comparing with the measured relative intensities and equiv- 
alent widths in our sample objects we can see that all of 
them are larger than the Starburst99 model predictions for 
a constant star formation law and, in some of them, like 
J1509 and J1540, even larger than the values predicted for 
an instantaneous burst. 

A great deal of work has been done in order to reconcile 
this disagreement between the observed luminosities of the 
WR bumps detected in the integrated spectra of galaxies 
and the model stellar atmospheres of WR stars and evolu- 
tionary synthesis models of clusters. Suggested hypotheses 
include the existence of a binary channel (Eldridge et al., 
2007) or the consideration of rotation in the stellar evolu- 
tionary tracks, which reduces the mass required to form WR 
stars (Meynet & Maeder, 2005). Nevertheless, some other 
reasons can be explored such as geometrical effects. If we 
consider the whole burst of star formation and we assume a 
continuous star formation law, we must compare the lumi- 
nosities of the WR bumps with the intensities of the Balmer 
lines and the continuum coming from the whole burst. This 
approach was already used by Perez-Montero & Diaz (2007) 
to match the observed fluxes of the WR features in Mrk 
209. On the contrary, if only the ionizing stellar population 
producing the WR stars is considered, only the burst region 
where the WR population concentrates must be analyzed. 
This could explain as well the lack of WR emission in the 
other three galaxies with an assumed continuous star forma- 
tion history. In the case of our sample, no aperture correction 
factors have been taken into account either for the Hq lu- 
minosities or the continuum contribution at the WR bumps 
wavelength, due to the compact nature of the objects. In- 
deed, the discrepancy factors between our Hq luminosities 
and those measured in the SDSS catalogue using a 3 arcsec 
fiber is not larger than 1.5 in any case and even smaller than 
1 in two of them (J1528 and J1729). Therefore, the aperture 
factor does not seem to be the cause of the abnormally high 
values found for the EW and the relative intensities of the 
WR features. 

The analysis of the relative number of WR to O stars 
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Figure 6. Relation between the WR equivalent widths and relative intensities of both the blue and the red bumps as a function of 
the age of the cluster, according to Starburst99 model predictions, at left for a metallicity of Z = 0.004 and at right for Z = 0.008. 
In all panels, black solid lines represent the instantaneous star formation model and red dashed lines the continuous star formation 
one. Crosses represent the corrected values derived for each galaxy and listed in Table 3 
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Figure 7. Left panels: time evolution of the ratio of WR. to O stars according to Starburst99 predictions. Right panel: same but 
for the ratio of WC to WN stars. In both cases, results for two mctallicities, Z= 0.004 and Z=0.008, are shown. Solid lines represent 
instantaneous star formation models and dashed line continuous star formation ones. Crosses represent the derived numbers for 
each galaxy taking into account the measured luminosities of the bumps as compared with those predicted by the Starburst99 
models at the predicted age. The number of O stars has been derived using the same procedure for L(Ha). 



Table 7. Luminosities of the Wolf-Rayct features measured for the sample of HII galaxies at 4650 A (blue bump, bb) and at 5808 A (red 
bump, rb), derived number of WR stars, ratios of WR to O stars and, for those galaxies with a detection of the red bump, of WC to WN 
stars. 



Object ID 


log L(bb) 
(erg ■ s-i) 


log L(rb) 
(erg • s-i) 


N(WR) 


WR/O 


WC/WN 


J0021 


39.84±0.09 




2500±600 


0.02±0.01 




J0032 


38.53±0.08 


37.79±0.16 


150±30 


0.06±0.01 


0.30±0.14 


J1509 


39.73±0.09 


38.84±0.11 


1750±400 


0.11±0.02 


0.19±0.08 


J 1540 


38.18±0.06 




60±15 


0.26±0.07 




J1616 


37.32±0.23 




9±6 


0.10±0.07 




J1624 


39.46±0.08 


38.92±0.12 


1070±200 


0.03±0.01 


0.64±0.22 


J 1729 


38.92±0.12 




300±100 


O.OliO.Ol 





Tailor-made models of HII galaxies 13 




0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 

t([SII]) (derived) 



0.9 1 1.1 1.2 1.3 1.4 1.5 1.6 

t([SIII]) (derived) 



0.9 1 1.1 1.2 1.3 1.4 1.5 1.6 

t([Sni]) (derived) 



1.1 1.2 1.3 1.4 1.5 

t([Oni]) (derived) 



Figure 8. Comparison between tiie electron temperatures (in units of IC* K) as derived from tiie observations and from the tailor- 
made photoionization models. Blaek cireles represent a comparison between a temperature directly derived in the observations, 
while white squares represent a temperature of another ion, assumed to be similar in the thermal structure. 



and of the number of WC to WN stars leads to very similar 
conclusions. In all cases the derived ratio of the number of 
WR to O stars is much higher than the asymptotic value 
reached with a continuous star formation, which is 0.022 for 
a metalhcity in between the two studied cases. For J1540, 
the obtained value, 0.44, is even much higher than the value 
expected for an instantaneous burst. As in the case of the 
luminosities of the WR features, we can reconcile the de- 
rived values enhancing the number of O stars on the basis 
that a non- negligible flux of the radiation from these stars 
is obscured by internal extinction. As in the previous case, 
probably geometrical effects can be determinant. Regarding 
the ratio between WC and WN stars, the derived values are 
lower in all cases than expected. This is due to the fact that 
the luminosities of the red bump predicted by the models 
are in much better agreement with the observed values than 
those of the blue bump and, hence, the relative number of 
WC stars is lower. 



4.3 Physical conditions and ionic abundances of 
the gas through tailor-made models 

The correct determination of ionic abundances in ionized 
gaseous nebulae using coUisional emission lines versus re- 
combination lines has become a matter of debate in the last 
years due to the disagreements predicted by photoioniza- 
tion models in the high metallicity regime {e.g. Chariot & 
Longhetti, 2001) or the discrepancies found in local nebulae 
{e.g. Peimbert et al., 2007 and references therein). 

Many times, these discrepancies are due to the sim- 
ple assumptions made about the ionization structure of the 
nebulae. Indeed, the only electron temperature measured in 
many cases is t([Oin]) because the emission line of [Oiii] at 
4363 A is the brightest auroral line. This temperature is as- 
sumed to be representative of the high excitation zone and 
the low excitation zone temperature is obtained under very 



simple assumptions. This can produce large deviations in 
high metallicity nebulae, where the low excitation ions have 
a high relative weight in the total abundance of a specific 
species (Perez-Montero & Diaz, 2005). One way to overcome 
this issue is to try to trace as precisely as possible the inner 
thermal structure of the gas. 

In Figure |B] we show a comparison between derived and 
modeled temperatures for the objects of our sample, for 
which at least four electron temperatures have been mea- 
sured. In the case of ions whose temperatures can not be 
derived due to the lack of observable lines in the optical, an- 
other line temperature, assumed to be representative of the 
zone where these iones are formed, has been used: t([Slll]) 
in the case of [Arm], t([Oni]) in the case of [Nelll] and [Felll] 
and t([Oll]) in the case of [Nil] in the objects where this was 
not available. As we can see, the agreement for t([Olll]) is ex- 
cellent, mostly due to the input conditions in the models to 
fit the auroral line of [Olll]. Regarding t([On]) and t([Sin]), 
the agreement is still good, although slightly worse that in 
the previous case. Finally, the largest deviations are found 
for the electron temperature of [Sll]. Except in the cases of 
J1455 and J1616, models predict much higher temperatures 
than observed. This is not surprising if we take into account 
that the emission of [Sll] is often found as well in the dif- 
fuse medium {e.g. Hunter & Gallagher, 1990) and therefore 
we are probably tracing in our observations a non-negligible 
portion of colder and less dense gas, which has not been 
considered in our models. This is compatible with the signi- 
ficatively lower radii found in our models in comparison with 
the values obtained from the spatial profile of Ha. Perhaps, 
this could be overcome if we considered in our models a low- 
density region in the outer parts of the nebula where lines 
from low excitation ions like [Sll] were produced. In the case 
of the three objects with a measurement of the electron tem- 
perature of [Nil] , the agreement is rather good in the case of 
J0021, while it is sensibly underestimated by the models in 
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Figure 9. Radial profile of the electron temperature, in violet solid line, and the relative abundance of different ionic species (from 
left to right and up to down, O, S, N, Ne, Ar and Fe) as calculated in the best photoionization tailor-made model of the galaxy 
J1616. We show as well the derived electron temperatures of oxygen and sulphur as black squares in the position of the average 
abundance of the corresponding ions. 
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Figure 10. Comparison between ionic abundances as derived from the observations using tlie direct method and the predictions 
of tailor-made photoionization models. 



the case of J1624 and J1729. This temperature is assumed 
to be equal to t([Oll]) in objects without a measurement of 
the auroral line of [Nil] at 5755 A. As we can see in the cor- 
responding panel, this assumption seems to work for most 
of the objects, except in the cases of J0032 and J1455. For 
these objects, we can not rule out some geometrical effect 
not considered in our photoionization models. 

For the rest of assumptions, as t([Arlll]) ~ t([Slll]) or 
t([Nelll]) ~ t([Felll]) ~ t([Olll]), we can see that all of them 
work fairly well. 

In Figure [9] we illustrate the good agreement found be- 
tween models and observations for the galaxy J1616. In each 
panel, in red, green and blue solid lines, the radial distribu- 
tion of the ionic fractions of the different elements as com- 
puted by the model (from left to right and from up to down: 
O, S, N, Ne, Ar and Fe) are plotted. The magenta thick solid 
lines show the radial distribution of the electron tempera- 
ture predicted by the same model. Finally, black squares 
show the electron temperatures obtained from the observa- 
tions. We have associated these points to the average zone 
where the corresponding ion is formed. In the cases of Ne2+ 
and Fe^+, t([Olll]) is shown; for Ar^+, t([Slll]) is plotted 
and in the cases in which no measurements of t([Nll]) are 
available, this is represented by t([Oll]). 

In Figure [To] we show the comparison between the rela- 
tive ionic abundances, derived using the corresponding mea- 
sured electron temperatures and the observed emission-line 
intensities, and the same abundances as predicted by our 
photoionization models. We can see that the agreement is 
somewhat better for high excitation ions, like 0^+, than for 
low excitation ones, like O"'' and S^. At the same time these 
latter ions present abundances spanning a wider range than 
the former. This might indicate that geometrical effects {e.g. 
density variations, matter bounded nebulae, etc..) affect the 
diagnostics of the gas in the low excitation region, hence a 
simple scheme can not be adopted for its analysis. In the 
case of our sample of HII galaxies, these differences are not 



appreciable due to their low metallicities and high excitation 
conditions, so the inner parts have a large weight on their 
average conditions. On the contrary, the analysis of disk HII 
regions and high metallicity objects are tied to a much larger 
uncertainty if these effects are not taken into account using 
as many indicators of the physical conditions of the gas as 
possible. 



4.4 The hardening of the radiation 

The stellar effective temperature has been traditionally the 
functional parameter most difficult to obtain. Vflchez & 
Pagel (1988) defined the so-called softness parameter, rj, 
which allows to estimate the hardening of the radiation field 
by means of the ratio of the ionic abundances of oxygen and 
sulphur, derived from the optical spectrum. 



V = 



o+/o^ 
s+/s^- 



(5) 



The same ratio, defined using the intensities of the cor- 
responding emission lines, can also be used to figure out the 
scale of equivalent effective temperature: 



V 



[OJJ]3727/[OJJ/]4959, 5007 
[5'//]6717, 6731/[Sm]9069, 9532 



(6) 



Both expressions agree in finding that HII galaxies 
present very high equivalent effective temperatures as com- 
pared to other HII regions of higher metallicity. Neverthe- 
less, there is a disagreement between the position of these 
objects in diagnostic diagrams and the scale of temperatures 
of the stellar atmospheres of O stars. HII galaxies usually 
lie in a region of the diagrams which corresponds to tem- 
peratures even higher than the maximum value of model at- 
mospheres, which is 50 kK. The cause has been sometimes 
attributed to some extra source of heating in HII galaxies 
(e.g. Stasihska & Schaerer, 1999). 
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Figure 11. At left, relation between the ratio of ionic abundances S+/S'^^ and O^/O^"*", whose quotient is the rj parameter. Black 
dots with error bars represent the modeled objects with the corresponding derived ratios, while blue squares are the values derived 
from the models. Red solid lines represent the fits to a sequence of photoionization models with the same effective temperature. 
The right panel represents the relation between the ratios of [SII]/[SIII] and [OII]/[OIII]. Symbols are the same as in left panel. 



The left panel of Figure [TT] shows the logarithmic ionic 
fractions 0+/0^+ vs. S+/S2+. The different lines of con- 
stant slope 1 in this diagram correspond to different values 
of the 7] parameter and hence of ionizing radiation temper- 
ature. The values derived from actual measurements of our 
observed objects are plotted as black solid circles while blue 
squares correspond to the values derived from the computed 
models. Red solid lines represent polynomial fits to different 
photoionization models computed using as ionizing source 
WM-Basic stellar atmospheres (Pauldrach et al., 2001) for 
different values of the effective temperature. These model 
sequences can help to establish the scale of the equivalent 
effective temperature. We can see that the models are not 
able to reproduce the high values of log(S^/S^"'"), and in 
some cases log(0'''/0^^) found in our objects that, in fact, 
would correspond to extremely high values of equivalent ef- 
fective temperatures. 

This disagreement disappears when the emission lines 
instead of the actual ionic ratios are considered. This is due 
to the fact that emission-line intensities are well fitted in all 
the tailor-made models while this is not the case for electron 
temperatures (see Figure[8]), specially for T([SII]) which is, 
in general overestimated by the models leading to an un- 
derestimate of the S^/H'^ abundance ratio. Again, the dis- 
agreement could be caused by the geometry of the diffuse 
gas in these HII regions, which is not adequately modeled, 
thus avoiding the need for an extra source of heating to ex- 
plain the apparently high stellar effective temperatures in 
these objects. At the same time, this probably means that 
the calibration of the effective temperature is not unique 
but in fact sensitive to the object class involved, through 
the different geometries of the ionized gas. 



5 SUMMARY AND CONCLUSIONS 

We have studied the underlying and ionizing stellar popu- 
lations and have modeled the properties of the emitting gas 



in a sample of HII galaxies. This sample was observed (H06 
and HOB) using double-arm spectrographs, which allow the 
simultaneous analysis of the entire spectral range from [Oil] 
3727 A to [SIII] 9532 A providing the derivation of the elec- 
tronic temperatures of, at least, four ions in all the cases: 
0"*", O^"*", S"*", and S^"*" and hence a precise determination of 
the ionization structure of the gas. 

The fitting of the blue part of the continuum was made 
using the program STARLIGHT and the SEDs computed 
from the Starburst99 stellar population synthesis models. 
This approach was used to study in a consistent manner the 
differences between the stellar populations in our sample of 
galaxies. The resulting distribution of ages and masses from 
this fitting revealed the presence of various stellar popula- 
tions. An underlying stellar population of about several Gyr 
dominates the stellar mass, while the blue to visual fight is in 
all cases dominated by young stellar populations with ages 
younger than 10 Myr. The presence of the older population 
can affect the reddening determination but its contribution 
to the measured EW(H/3) is only about 10% for most of the 
sample objects, being around 25% in the case of J0032 and 
almost zero for J1616 and J1729. In fact, regarding the sub- 
traction of Balmer absorption features from the underlying 
stellar populations, we did not find that the spectral fitting 
technique constitute a substantial improvement over the si- 
multaneous fitting of absorption and emission components 
as performed in previous works (HOG and H08). 

We analysed the Wolf-Rayet clusters detected in seven 
of the ten studied galaxies through the WR broad features 
seen in their spectra. The fitting of the stellar population 
and the modeling of the ionized gas allows to correct them 
for the contribution of the underlying population and the 
nebular emission to the continuum and for the dust absorp- 
tion factor in the case of the emission of the Balmer lines. 
The comparison between the corrected fiuxes and equiva- 
lent widths of the WR blue bump with the predictions from 
Starburst99 models at the corresponding metallicity leads to 
the estimate of the number of WR stars which ranges from 
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9±6 to 2500±600 depending on the galaxy. In three of them, 
also the red bump was measured allowing to estimate the 
WC to WN number ratio, Nevertheless, the fluxes of both 
the blue and red bumps relative to H/3 are larger than the 
values predicted by the models in almost all the cases. 

A continuous star formation history of the stellar ion- 
izing population has been assumed in all the objects of the 
sample. This is supported by the simultaneous presence of 
very massive O and WR stars, along with the relatively low 
measured values of EW(H/3), once corrected for underlying 
contribution and dust absorption. The ages of the star for- 
mation bursts, the number of ionizing photons as derived 
from the Ha fluxes and the intensities of the most promi- 
nent emission lines were used to constrain photoionization 
models for each galaxy. The Starburst99 libraries were used 
to calculate the SED of the ionizing cluster in order to be 
consistent with the previous fitting of the stellar popula- 
tions. Some dust is assumed to be present associated with 
the gas. The amount of dust in each model was derived in 
order to reproduce the observed EW(H/3). These derived 
values agree well with the main observed quantities in each 
object, with the exception of the sizes of the ionized regions 
which result 5-10 times smaller in the models as compared 
to measurements from Ha in the spatial profiles of the slits. 

The observed electron temperatures and ionic abun- 
dances are, in most cases, well reproduced by the models 
which implies that the direct method using coUisional emis- 
sion lines provides robust answers if the thermal structure of 
the nebula is well traced. Therefore, our models were able to 
reproduce the ionization structure of the nebula without ap- 
pealing to temperature fluctuations or any other geometrical 
effects that can affect the abundances derived from collision- 
ally excited emission lines. This result is consistent with the 
temperature measured from the Balmer jump in some of the 
objects in our sample (H06) . Only the electron temperature 
and the ionic abundances of S"^/H+ is not well reproduced 
by our models, which could be consequence of the presence 
of diffuse gas in these galaxies, not taken into account in 
the models. This would also be consistent with the small 
sizes predicted by the models as compared with Ha spatial 
profiles in the slits. 

Finally, regarding the hardness of the radiation field, 
model predictions agree with observations when the soft- 
ness parameter rj, which parametrizes the temperature of 
the ionizing radiation, is expressed in terms of emission line 
intensity ratios. However, there is a disagreement when ionic 
abundance ratios are used instead. This is probably caused 
by the overestimate of the electron temperature of by the 
models and the corresponding underestimate of the S'''/S^^ 
abundance ratio rather than by existence of an additional 
heating source in HH galaxies. 
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APPENDIX A: TOTAL ABUNDANCES AND 
IONIZATION CORRECTION FACTORS 

The calculation of tailor-made photoionization models for 
each galaxy with a reliable inner ionization structure based 
on the measurement of various electron temperatures al- 
lows a complete description of the total abundances of the 
main elements present in the gas. With this information it is 
possible to compare the ionization fractions obtained from 
the models with the most widely used ionization correction 
factors (ICFs), necessary to derive the total abundances of 
certain elements for which there are no lines in the optical 
range for all the main ionization stages. Generally, for an 
ionic species X'^\ whose abundance has been derived from 
forbidden coUisional lines, we define the ICF for that ion in 
the following manner: 



(Al) 



In Table IATI we list the total abundances and the ICFs 
obtained from the models for N, S, Ne, Ar and Fe. 

The oxygen ionic abundance ratios, O"*"/!!"^ and 
0^+/H+, usually derived from the [On] AA 3727,29 A and 
[Qui] AA4959, 5007 A lines are related to the total abun- 
dance of oxygen assuming that: 

O _ 0+ + 
H ^ H+ 



which is based on the assumption that the relative abun- 
dances of neutral hydrogen and oxygen are similar. We have 
checked in our models that this assumption is valid within a 
1% of uncertainty. As the bright emission lines of oxygen are 
well reproduced by our models, as also is the temperature 
sensitive line, the agreement between the total oxygen abun- 
dances derived from observations and models agree fairly 
well with deviations, in most cases, no larger than the re- 
ported observational errors. 

The total abundance of nitrogen is obtained from the 
ionic abundance of N+, derived from the intensities of the 
AA6548, 6584 A lines. Then, a similar ionization structure 
for nitrogen and oxygen is assumed, according to the follow- 
ing expression: 

iV+ 

0+ 



O 



(A3) 



which, combined with equation (6), leads to the following 
ICF for nitrogen: 



ICF(N^ 



1 + 



91 



(A4) 



In the top left panel of Figure lAll we plot these ICFs 
as a function of the O^/O^^ ratio as compared to the best 
linear fit to the models. The slope of the fit is slightly higher 
(1.33), than the value derived using the classical assumption. 
This is a consequence of the fact that the ionic fraction of 
N""" is lower than that of in these type of objects. 

Regarding sulphur, we have derived S^/H^ abundances 
from the fluxes of the [Sll] emission lines at AA 6717, 6731 A 
and T([Sll]) and S'^'^/H^ abundances from the fluxes of 
the near-IR [Slll] AA 9069,9532 A lines and T([Slll]). The to- 
tal sulphur abundance has been calculated using an ICF 
for S'''-|-S^''' according to Barker's (1980) formula, which is 
based on the photoionization models of Stasihska (1978). We 
have written this formula as a function of the ratio 0^"^/0 
instead of 0+/0 in order to reduce the propagated error for 
this quantity. 



ICF{S+ + 



6' 



0+ + 



-1/cy 



(A5) 



where a — 2.5 provides the best fit to the scarce observa- 
tional data on S"^^ abundances (Perez-Montero et al. 2006, 
PM06 in the plot). Nevertheless, we can see in Figure [ATI 
that a value of a = 3.3 provides a better fit to the results 
obtained from the photoionization models. 

Neon is only visible in the spectra by means of the 
[Nelll] emission line at A3869 A. The total abundance of neon 
has been calculated using the following expression for the 
ICF (Perez-Montero et al., 2007, PM07 in the plots): 



ICF{Ne^ 



0.7534-0.142- 



0+ + 02+ 



Q+ I q2+ 

+0-171 02+ (A6) 



(A2) 



This formula takes into account the overestimate of 
Ne/H in objects with low excitation, where the charge trans- 
fer between O^^ and H" becomes important (Izotov et al., 
2004, 104 in the plot). Nevertheless, the total abundances ob- 
tained in our models are slightly higher than those derived 
with the use of this ICF. 

The only available emission lines of argon in the op- 
tical spectra of ionized regions are those of [Arm] at 
7136 A and [Ariv] at 4740 A. It is usually assumed that 
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Table Al. Ionization correction factors and total abundances as derived from photoionization models of the main elements with emission 
lines in the optical spectral range of the studied HII galaxies. 





J0021 


J0032 


J1455 


J1509 


J1528 


J1540 


J1616 


J1624 


J1657 


J1729 


12+log(0/H) 


7.99 


8.01 


7.94 


8.18 


8.18 


8.15 


8.00 


8.04 


8.01 


8.12 


ICF(N+) 


4.17 


4.79 


6.61 


5.13 


3.80 


3.47 


9.78 


5.13 


3.39 


5.37 


12+log(N/H) 


7.13 


6.78 


6.76 


7.01 


7.00 


6.90 


6.62 


6.71 


6.76 


7.16 


ICF(S++S2+) 


1.11 


1.14 


1.28 


1.12 


1.10 


1.04 


1.37 


1.16 


1.07 


1.19 


12+log(S/H) 


6.29 


6.40 


6.23 


6.57 


6.41 


6.54 


6.33 


6.30 


6.34 


6.41 


ICF(Ne2+) 


1.20 


1.17 


1.10 


1.20 


1.23 


1.29 


1.07 


1.15 


1.26 


1.15 


12+log(Ne/H) 


7.24 


7.24 


7.16 


7.36 


7.43 


7.22 


7.17 


7.22 


7.19 


7.37 


ICF(Ar2+) 


1.15 


1.17 


1.23 


1.15 


1.15 


1.10 


1.29 


1.17 


1.15 


1.17 


ICF(Ar2++Ar3+) 


1.05 


1.06 


1.03 


1.04 


1.07 


1.06 


1.02 


1.04 


1.08 


1.04 


12+log(Ar/H) 


5.58 


5.69 


5.58 


5.87 


5.78 


5.82 


5.64 


5.68 


5.57 


5.79 


ICF(Fe2+) 


2.69 


2.88 


3.89 


2.95 


2.57 


2.14 


5.13 


3.09 


2.40 


3.24 


12+log(Fe/H) 


5.82 


5.80 


5.86 


5.92 


5.99 


5.73 


6.15 


5.77 


5.79 


5.88 




0'*/(0*+0'*) S'*/(S*+S'*) 0'*/(0*+0'*) s'*/(S*+s'*) 



Figure Al. Comparison between different ionization correction factors (ICFs) derived for different elements from the computed 
photoionization models (black circles) and some of the most used relations. 



T([Arlll])«T([Sin]) (Garnett, 1992). The total abundance 
of Ar can be calculated using the ICF(Ar^+) derived from 
photo-ionization models by Perez-Montero et al. (2007, 
PM07 in the plots), in the case that only the emission line 
of [Arm] be detected 

ICF{Ar'^+) = 0.749 + 0.507 • (l - ^f^, t,.,^ ) + 



0+ + 02+ 



+ 0.0604 ■ 1 - 



0+ + 02+ 



(A7) 



In the same work, another ICF is proposed for the case in 
which [Ariv] is also measured: 



ICF{Ar+ + Ar+) = 0.928 + 0.364 • ( 1 



O 



2+ 



0+ + 02+ 



+ 



+ 0.006 • 1 



i2+ 



o 

0+ + 02+ 



(A8) 



These two ICFs can be expressed as well as a function 
of the ionic abundances of sulphur, and S^^, when only 
red spectroscopic observations are available: 



ICF{Ar'^+) = 0.596 + 0.967 • ( 1 



^2+ 



S+ + S2+ 



+ 0.077 • 1 



S+ + S2^ 



ICF{Ar'^+ + Ar^+) = 0.870 + 0.695 • ( ( 1 



s+ + 



(A9) 



+ 



+ 0.0086 ■ 1 - 
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In the bottom panels of Figure lAll we can see that the 
agreement between these expressions and the ICFs found in 
our models is good. 

Finally, iron abundances in the gas-phase can be cal- 
culated by correcting the ionic abundances of Fe^"*" , derived 
from the intensity of the emission line of [Felll] at 4658 A, us- 
ing the following expression from Rodriguez & Rubin (2004, 
RR04 in the plot): 



with a value of a = 0.09. Nevertheless, the fit to the ICFs 
found in our models gives a best value for a — 0.24, which 
leads to lower total abundances of iron. 




(All) 



